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Abstract 

The metal-free (Pop. Ill) and extremely metal-poor (EMP) stars of low- and 
intermediate-masses experience mixing of hydrogen into the helium convection dur- 
ing the early TP-AGB phase, differently from the meal-rich stars. We study the 
nucleosynthesis in the helium convective zone with 13 C formed from mixed protons 
as neutron source by using a nuclear network from H through S. In the absence or 
scarcity of the pristine metals, the neutron-recycling reactions, 12 C(ri,7) 13 C(a,r2,) 16 
and also 16 O(n,7) 17 O(a,n) 20 Ne promote the synthesis of O and light elements, in- 
cluding their neutron-rich isotopes and the odd atomic number elements. Based on 
the results, we demonstrate that the peculiar abundance patterns of C through Al ob- 
served for the three most iron-deficient, carbon-rich stars can be reproduced in terms 
of the nucleosynthesis in Pop. Ill, AGB stars in the different mass range. We argue 
that these three stars were born as the low-mass members of Pop. Ill binaries and 
later subject to the surface pollution by the mass transfer in the binary systems. It is 
also shown that the AGB nucleosynthesis with hydrogen mixing explains the abun- 
dances of C, O, Na, Mg and Al observed for most of carbon-enhanced EMP (CEMP) 
stars, including all CEMP-s stars with s-process elements. In addition the present 
results are used to single out other nucleosynthetic signatures of early generations of 
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stars. 
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1. Introduction 

During the past decade, the studies of extremely metal-poor (EMP) stars in the Galactic 
halo with high resolution spectroscopy have revealed peculiar abundance patterns, different from 
those known from previous observations and predicted from theoretical calculations. Striking 
examples are the two stars of giant HEO 107-5240 (Christlieb et al., 2002) and sub-giant HE1327- 
2326 (Frebel et al., 2005; Aoki et al., 2006), which have the smallest iron abundances among 
the known EMP stars, i.e., of [Fe/H] = —5.3 and —5.4, respectively. The both stars share the 
characteristic feature among the light element abundances; apart from the huge enrichment of 
carbon ([C/Fe] =4.0 and 4.26), the light elements such as nitrogen, oxygen, sodium, magnesium, 
and aluminum are reported to be largely enhanced (Christlieb et al., 2002, 2004; Bessell et al., 
2004; Frebel et al., 2005, 2006; Aoki et al., 2006). Some differences are also discerned between 
these two stars; the enhancements of nitrogen, oxygen, and Na through Al are much larger 
for HE1327-2326 than for HE0107-5240 (by 2.2 dex, 0.5 dex, and 1.2-1.5 dex, respectively). 
Moreover, the former displays the enrichment of an s-process element, strontium ([Sr/Fe] = 
1.1 ~ 1.3), while the latter does not ([Sr/Fe] < —0.5). Recently, another giant HE0557-4840 is 
identified at the iron abundance of [Fe/H] = —4.75, much smaller iron abundance than other 
EMP stars of [Fe/H] > —4 (Norris et al., 2007), which also has a large carbon enrichment but 
smaller than the above two stars ([C/Fe] = 1.65 — 2.08). 

In the theoretical point of view, new or improved models of stellar evolution and nucle- 
osynthesis are required to explain such extraordinary abundance patterns observed, including 
both the similarities and differences. In this paper, we call the two stars of [Fe/H] < —5 and 
the other star of [Fe/H] < —4 as hyper metal-poor (HMP) and ultra metal-poor (UMP) stars, 
respectively, as proposed by Beers & Christlieb (2005) to distinguish them from the other 
EMP stars of [Fe/H] > —4. Several theoretical approaches have been made to understand the 
abundance patterns observed for these two HMP stars. A straightforward interpretation of 
the observed small iron abundances is that these stars are the second-generation stars, born 
from gas polluted by the ejecta of first generation supernova. Umeda & Nomoto (2003) and 
Iwamoto et al. (2005) proposed a peculiar supernova model of metal-free 25 M Q star that gives 

* JSPS Research Fellow. Present Address : Astronomical Data Archives Center, National Astronomical 
Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 

t Institut d'Astronomie et dAstrophysique, C.P.226, Universite Libre de Bruxellcs, B-1050 Brussels, Belgium 

* Astrophysics Group, EPSAM, Keele University, Keele, Staffordshire ST5 5BG, UK 



2 



an ejecta of unusual abundances with little iron and large amount of CNO elements, similar 
to those observed for HE-0107-5240 and HE1327-2326, respectively (Tominaga, 2009, see also). 
Limongi et al. (2003) presented a scenario involving two metal-free supernovae, one with low 
mass (~ 15 M ) and the other with high masses (~ 35 M ) for HE0107-5240. Meynet et al. 
(2006) argued a combination of stellar wind and supernova ejecta of rotating massive star for 
HE1327-2326. For all these scenarios adopt the crude supernova ejecta, however, the oxygen 
abundance may be a critical observation in discriminating possible scenarios as pointed out by 
Bonifacio et al. (2003). Actually, Bessell et al. (2004) and Frebel et al. (2006) derive rather 
small oxygen enrichment in disagreement with these supernova scenarios. 

The EMP stars, currently observed in the Galactic halo, are the low-mass stars and have 
survived their long lives to date. We should take into account the modifications of their surface 
abundances through not only the internal but also the external processes after their birth. As 
discussed in the previous paper (Suda et al., 2004), the origin of iron-group elements with 
such tiny abundances can readily be explained in terms of the accretion of the interstellar gas, 
enriched with metals ejected from the supernovae of later generations, in the clouds where they 
were born. In addition, it is also argued that the enhancements of carbon and light elements are 
consequent upon the surface pollution by mass transfer from an AGB companion in the binary 
systems of metal-free stars. The latter possibility attracts our interest because it suggests that 
these stars may belong to the first generation stars, which are called as Population III (Pop III). 
The stars that have formed just after the Big Bang can provide information on the formation of 
first-collapsed objects and on the chemical and dynamical evolution in the very early universe. 

In order to discuss the possibility that the HMP and UMP stars are low-mass survivors of 
Pop III stars, we have to develop new aspect of evolution of extremely metal-poor stars. For the 
metal-free and extremely metal-poor (EMP; [Fe/H] —2.5) stars, it is proved that the helium- 
flash convection extends into the hydrogen-rich envelope in the stars of mass M < 3.0 M , 
different from metal-rich populations (Fujimoto et al., 1990, 2000). What causes this difference 
is a lower entropy of the hydrogen-burning shell in the stars of smaller CNO abundances, 
and hence, a smaller entropy barrier between the helium-zone and the tail of hydrogen-rich 
envelope. The engulfed hydrogen is carried down by convection and captured by 12 C in the 
middle of helium convective zone to create 13 C through the /3-decay of 13 N. Then 13 C is 
carried further inward to emit the neutrons through the 13 C(a,ri) 16 reaction, which initiates 
the neutron capture reactions and the s-process nucleosynthesis. The nuclear products in the 
helium convective zone will later be brought up to the surface, dredged up by the envelope 
convection that deepens in mass as the burning shells expand absorbing the energy released 
during the shell flashes (Hollowell et al., 1990; Fujimoto et al., 2000; Suda et al., 2004; Iwamoto 
et al, 2004). 

The nucleosynthesis triggered by the 13 C-burning in the helium-flash convection was orig- 
inally proposed by Sanders (1967), relying upon an early computation of helium shell flashes 
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by Schwarzschild & Harm (1967), which met with hydrogen mixing into the helium convection 
for 1M stars of Population II composition. Ulrich (1973) elaborated this line of thought to 
propose a plum mixing model and argued the repeated irradiation of neutrons during the recur- 
rence of helium shell flashes (see also Scalo & Ulrich, 1973). But it has been later demonstrated 
that the entropy barrier becomes too high and prevents the helium flash convection from reach- 
ing the tail of hydrogen-rich envelope if the radiation pressure is properly taken into account 
in the Pop. I and II stars (Iben, 1976; Fujimoto, 1977). Iben & Renzini (1982) then argued 
that hydrogen could be carried down to form 13 C by the semi-convection which developed in 
the carbon-rich helium zone during the third dredge-up ; they propose that the latter acts as 
neutron source as it is incorporated into the helium convection during the next shell flash, but 
later, Straniero et al. (1995) found that 13 C burned during the inter-pulse phase. On the other 
hand, it is shown that for very thin envelope, the entropy of envelope decreases to allow the 
helium flash convection to extend through hydrogen-rich layer (Fujimoto, 1977). Finally, it is 
revealed that in the extremely metal-poor stars, the entropy of hydrogen-rich envelope can be 
small enough to mix hydrogen into the helium flash convection (Fujimoto et al., 1990, 1995, 
2000). 

The purpose of the present work is to explore the characteristics of nucleosynthesis in 
the helium-flash convective zone, triggered by the hydrogen injection in the extremely metal- 
poor stars of low and intermediate masses. It has been suggested that the neutron-recycling 
reactions, 12 C(n, 7) 13 C(a, n) 16 0, play a key role in producing oxygen (Gallino et al., 1988; 
Jorissen & Arnould, 1989) since neutrons, emitted in the helium-flash convection, may be 
absorbed by 12 C because of the large abundance. A definite conclusion could not be reached in 
those days, however, because of large uncertainties of neutron capture cross sections of 12 C and 
16 0. Recently, the neutron capture rates for these nuclei have been improved (Ohsaki et al., 
1994; Igashira et al., 1995, and also see Bao et al. 2000), and this enables us to draw a more 
quantitative and sound picture. Under the initial scarcity of neutron poisons in the extremely 
metal-deficient stars, the neutron recycling reactions are expected to have critical importance, 
and may further lead to the generation of heavier elements through the sequences of neutron and 
a-captures. In the present paper, we focus on the progress of neutron recycling reactions and 
the resultant formation of light elements from oxygen through sulfur including their isotopes 
as well as odd number nuclei of sodium, aluminum and phosphorus. The synthesis of heavy 
elements beyond sulfur and of the s-process elements will be discussed in a subsequent paper 
(Nishimura et al., in preparation). Since the nuclear products, thus synthesized in the helium 
convection, are eventually carried out to the surface, the present results may have relevance 
to the observed abundances for the extremely metal-poor stars. In particular, the comparison 
with the abundance patterns observed from the two most metal-deficient stars, HE 0107-5240 
and HE 1327-2326 may give an insight into their origins. 

The convective 13 C-burning model, studied in this paper, is different from the so-called 
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radiative C-burning model. The former is founded on the evolution of EMP stars within the 
current standard framework of stellar evolution. In contrast, the latter resorts to the alleged 
formation of 13 C-pocket in the radiative helium zone during the third dredge-up; hydrogen 
is assumed to be injected from the bottom of surface convection into the top of helium zone 
where a thin 13 C-rich layer is formed via the proton-capture of 12 C. The formation of 13 C- 
pocket is presumably due to the convective overshooting, though it has not yet established. If 
it is properly formed, 13 C burns to emit neutrons and trigger the s-process nucleosynthesis in 
the radiative zone at relatively low temperatures before the next helium flash ignites (see the 
review by Busso et al., 1999). It has been attempted to apply the radiative 13 C-burning model 
to the s-process nucleosynthesis in a zero-metal AGB stars (Goriely & Siess, 2001) and in very 
metal-poor AGB stars (Z = 5 x 10" 5 , Straniero et al., 2004). 

We describe the methods and assumptions of the present work in the next section. 
Results of the nucleosynthesis in the helium flash convective zone, computed with the adopted 
nuclear reaction network, are presented and examined in Section 3. In Section 4, the character- 
istics of their yields are examined and compared with the observations for HMP / UMP stars 
and EMP stars on the production of carbon and oxygen and also on the relative abundances 
of sodium, magnesium, and aluminum. Then, we inquire into the nature of currently known, 
three EMP/UMP stars. The conclusions and discussion follow. 

2. Methods and Approximations 

A distinctive feature in the evolution of extremely metal-poor stars ([Fe/H] < —2.5) is 
that the helium convection penetrates into the hydrogen-rich layer and mixes hydrogen down 
into the convective zone (Fujimoto et al., 2000). This occurs for the stars of mass M<3M . The 
progress of nucleosynthesis and material mixing during the helium shell flashes is schematically 
illustrated in Figure 1. As the TP- AGB phase starts, the helium shell flashes recur and grow in 
strength to drive convection more outward, until it reaches the hydrogen-rich layer (Fujimoto 
et al., 1990, 2000). Hydrogen mixed into the helium-flash convection is captured by 12 C to 
form 13 C via the reactions 12 C(p, 7) 13 N(e + z/) 13 C in the middle of convective zone. It diffuses 
further into deeper interior and releases neutron via 13 C(a, n) 16 reaction to promote the 
nucleosynthesis via neutron capture reactions. These nuclear products are stored in the helium 
zone during the recurrence of shell flashes. Finally, when the helium flash grows strong to erode 
hydrogen at a sufficiently large rate, the burning of mixed hydrogen ignites a hydrogen-shell 
flash and splits the helium convective zone into two in the middle. This brings about the helium 
flash-driven deep mixing (He-FDDM) that conveys the nuclear products, stored in the helium 
zone and processed by hydrogen burning, to the surface (Hollowell et al., 1990; Iwamoto et al., 
2004). After He-FDDM enriches the envelope with CN elements above [CNO/H] ~ -2.5, the 
hydrogen mixing can no longer take place because of a large entropy barrier between the helium 
flash-convective zone and the bottom of hydrogen-rich envelope. If the star is more massive than 
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say, M ~ 1.5 M Q , the third dredge-up follows to bring up to the surface the nucleosynthetic 
yields in the helium flash-convective zone, including freshly synthesized 12 C. The s-process 
elements are also dredged up to the surface, as will be discussed in the subsequent paper. The 
detailed description about the model is given in the previous papers (Fujimoto et al., 2000; 
Suda et al., 2004; Komiya et al., 2007). 

We investigate the progress of nucleosynthesis in the helium-flash convection considering 
neutron production through the hydrogen mixing. The variations of the strength of shell flashes 
are provided by the simulations of hydrogen mixing performed with stellar evolution codes for 
various initial masses and pristine metallicity (see Suda & Fujimoto, 2008, and the references 
therein and Campbell & Latanzio 2008). Furthermore, the efficiency of elemental mixing is 
subject to uncertainties due to the treatment of convection and the entrainment of matter by 
convective elements. Since the purpose of this work is to explore the general characteristics of 
nucleosynthesis triggered by hydrogen mixing, we may well treat the amount of mixed hydrogen 
and the duration of mixing event as parameters. In our computations, we deal with the mixing 
of 13 C, instead of hydrogen, assuming that all the engulfed protons are captured by 12 C in the 
middle of convective zone. This is a good approximation as far as the number ratio of engulfed 
protons to 12 C is sufficiently small. The mole abundances of CN elements resultant from proton 
captures are estimated at 

F 13 /AX mix = 1 - a( AX mix /Y 12 ) (1) 
F 14 /AX mix = a(AX mix /F 12 )/2 (2) 

where AX m ; x is the amount of mixed hydrogen per unit mass and a is the ratio of proton 
capture rates between 13 C and 12 C [a = (cru)p 7 ( 13 C)/ (av ) P7 ( 12 C) ~ 3 ~ 4]. Accordingly we may 
well ignore the fraction of mixed hydrogen that is captured by 13 C to form nitrogen, as long as 
AX m i x /Y"i2 <C 1, i.e., the amount of mixed proton is much smaller than the amount of 12 C in the 
helium convection. In addition, 14 N, if mixed, has little effect on our results essentially since it 
promptly reacts with neutron via the 14 N(n,p) 14 C reaction and proton is absorbed mostly by 
12 C to form 13 C. 

The development of shell flashes is calculated with use of the analytic solution of burning 
shell by Sugimoto & Fujimoto (1978), which enables the coverage of wide parameter spaces. We 
further adopt the one- zone approximation, as formulated by Fujimoto (1982, see also Fujimoto 
et al. 1999) under the assumption of the uniform distributions of abundances in the helium 
convective zone and with use of the averaged nuclear reaction rates over the convective zone 
estimated under the approximation of adiabatic gradients of density and temperature. Figure 2 
illustrates the time variations in the temperature and density at the bottom of helium convective 
zone during the shell flashes. The model parameters are summarized in Table 1. The mass, 
M c , and radius, r c , of the core interior to the helium burning shell and the proper pressure P* 
at the bottom of helium burning shell are adopted from an evolutionary computations of EMP 
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stars and metal-free Pop. Ill stars when they experience the hydrogen mixing for the first time. 
Model 1 and 3 are from the evolutionary calculation of Pop. Ill models of 1.5M and 2.0M Q by 
Suda & Fujimoto (2008). Model 4 is from the evolutionary calculation of 2M and [Fe/H] = —3 
by Fujimoto et al. (2000). In addition, we devise Models 2 and 5 by slightly modifying Models 
3 and 4, respectively, in order to cover the parameter range and in particular, the latter to see 
the effect of large temperature. 

In the present study, we consider the amount, AY"i3 jm i x , of mixed 13 C in mole per unit 
mass as parameter, setting AY\ 3 iTT1 i Y = AJ^, and add it at a constant rate for an interval 
At mix starting from time t pea k when the burning rate of 3a reaction reaches the maximum. We 
may express the amount of mixed 13 C by the mass ratio, A( 13 C/ 12 C), to the 12 C abundance at 
the start of mixing, i.e., A( 13 C/ 12 C) = AX 13 , mix /X 12 (t pcak ) = (13/12) AF 13 , m i x /Fi 2 (Wk)- The 
amount of mixed 13 C corresponds to the amount of mixed hydrogen- rich matter as 



AM H , mix = 4.6xlO- 6 Hc - conv 1 w ^ 12 ' Hc - conv M , (3) 




where AMne-conv is the mass in the helium convective zone. Accordingly only a small part of 
the mass in the tail of hydrogen profile suffices to account for the amount of mixed hydrogen- 
rich matter. For A( 13 C/ 12 C) <C 1, this amount of mixed hydrogen is significantly smaller than 
reported for the stellar evolution models which develops He-FDDM (Hollowell et al., 1990; 
Fujimoto et al., 2000; Iwamoto et al., 2004; Straniero et al., 2004), and liable to be overlooked 
in the evolutionary computations. As for the ignition of hydrogen shell-flash and splitting of 
helium flash convection, Sweigart (1974) argues that it is related to the energy deposition rate 
due to the mixed hydrogen, and hence, to the mixing rate rather than to the amount of mixed 
hydrogen, and estimate dM^ m i x / dt ^3x 1O~ 14 M s _1 for the shell flash with the peak energy 
generation rate due to the helium flash L p = 1.2 x 1O 7 L (and the core mass M c = 0.62 M ). 
This is transferred to the lower limit to the duration of mixing epoch as 

Aimix > i.6 x 10 » fW) (AMh^\ (4) 

\10 7 L Q J V °- 01 / V 0.03 M J ' K J 

which is sufficiently smaller than the timescales of temperature decrease after the peak (see 

Figure 2 and Table 1). 

Along the temperature and density histories of the models, we follow the progress of 

nucleosynthesis by adding a specific amount of 13 C to the helium convective zone in a specific 

interval after the peak of helium shell-flashes. The initial abundances of helium, 12 C and 16 

in the helium convective zone are also taken from the models of metal-free stars. As compared 

with them, the other metals have much smaller abundances in the helium zone of EMP stars 

and their effects will be discussed later based on the results. The nuclear reaction network 

involves 59 isotopes of the 16 lightest species from the neutron and 1 H through 35 S, which is an 

extension of the version by Aikawa et al. (2001, see also Aikawa et al. 2004). Reaction rates are 

taken from NACRE (Angulo et al., 1999) and Caughlan & Fowler (1988) for charged-particle 
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reactions and Bao et al. (2000) for neutron capture reactions. As for the neutron reactions of 
17 0, the cross section of 17 0(n, a) 14 C is set at 10 mb, given by Koehler & Graff (1991) in the 
temperature range relevant to our problems. The cross section of 17 0(ri,7) reactions, which 
is missing out of the compilation by Bao et al. (2000), is set equal to that of 16 0. The (n, a) 
reaction of 33 S and the cross sections of (n, 7) reactions of 32 P are taken from Rauscher & 
Thielemann (2000). 

Figure 3 outlines the main paths of nucleosynthesis from 12 C through 35 S via the a and 
neutron captures. Under the presence of neutron source, the (a,n) reactions coexist with the 
(a, 7) reactions. In particular, the j4 Z(n,7) j4+1 Z(a,n)" 4+4 (Z + 2) reaction may compete with 
the A Z(a,7) j4+4 (Z + 2) reaction to form a-nuclei of mass number A + 4 from a-nuclei of mass 
number A. This competition may be expressed in terms of the ratio, Aa+4, of production rates 
of A+i {Z + 2) nuclei at the moment t, defined as: 

A m = Y A+1 (t)Y a (t)p(t)N A (av) an ( A ^Z) 

M+4[t) Y A (t)Y a (t)p(t)N A (av) a ,( A Z) ' l ° j 
which represents the dependence of the production of the isotopes on the reaction paths. Here, 
Y A (t), p(t) and N A (av) a b are the number abundance in moles, the density and the rate of (a,b) 
reaction, respectively. To investigate the influence of neutrons, Y A+ i(t) is assumed to be created 
only through A Z(n,'y) A+1 Z reaction without the destruction by the (a,n) reaction. We denote 
the production rate ratio under this assumption by A*. For simplicity, we assume 

Y n (t') = (t'< t mix ) 
Y A {t')/Y A {t)~\ (t'>t mix ) (6) 
(crv) rn ( A Z) ~a ni ( A Z)v n 

where t m { x means the moment when the mixing starts. The first and second equations assume 
that there are neither neutrons nor A+1 Z before the mixing and little changes in most of a- 
elements after the mixing, respectively. In the third equation, the cross section and the thermal 
velocity, v n , of neutrons are assumed to be constant instead of (av) ny . Then we can get the 
simple equation to estimate the effect of neutrons, 
\* U\ n-U\„ tA^ N A {av) an { A+l Z) 

AW*) = <t)a ni ( Z) Na{(TvU{Az) , (7) 

where the exposure r(t) is defined as 

r(t)= fY n {t')p{t')N A v n dt'. (8) 
Jo 

The above assumptions will not strongly affect our discussion because in almost all the cases 
A^ +4 (t) results in such large values as vary by several orders of magnitude (see below). 

Figure 4 shows the comparison of the (0,7) and (a,n) reaction rates for the light ele- 
ments. The A+1 Z(a,n) A+4 Z rates are much larger than the A Z{a^) A+i Z rates. Their differences 
amount to many orders of magnitude for carbon and oxygen. Thus, A^ +4 (t) ^> 1 so that the 
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(n, 7) and (a, n) reactions dominate the formation of a-nuclei even under small neutron ex- 
posures. It should be noted, however, that the actual ratio of the formation rates results 
much smaller than X* A+4 (t) if we take into account the destruction of A+1 Z nuclei by a-capture 
[\A+4(t) ~ (dr I 'dt)a ni ( A Z) / p(t)N a(cv) a7 ( A Z)]. For magnesium and beyond, the a-capture re- 
actions themselves may reduce to be too slow to occur during the helium shell flashes. For 
heavier nuclei, on the other hand, the cross sections of neutron capture grow larger in general. 
Accordingly, with under a small neutron exposure, the helium burning can drive the nucleosyn- 
thesis up to heavier nuclei than otherwise expected, via the combinations of neutron capture 
and (a,n) reactions. 

3. Results 

We compute the progress of nucleosynthesis in the helium convective zone during the 
helium shell flashes in the metal-free stars, born of the gas totally lacking metals, by adding 
13 C for the thermal history given by Models 1-5 in Figure 2. The progress of nucleosynthesis 
is exemplified in Figures 5, where the variations of the abundances are plotted for the major 
elements and isotopes against the elapsed time from the onset of shell flash (left panel) and 
against the neutron exposure (right panel) for Model 4 with A( 13 C/ 12 C) = 10~ 2 and At mix = 10 9 
s. As soon as the mixing starts, the 13 C(«,n) 16 reaction proceeds in much shorter timescale 
than 3a reaction and the neutron density soon reaches as large as N n ~ 10 11 cm -3 . Because of 
the lack of pristine metals as absorbers, neutrons are exclusively react with the most abundant 
nucleus 12 C and enter into the neutron-recycling reactions of 12 C(n, 7) 13 C(a, n) 16 (Gallino 
et al., 1988). Along with the run of neutron- recycling reactions, the second most abundant 
nucleus 16 increases in number, eventually to absorb a neutron and to cause the leak from the 
recycling circuit. 

The formation of 17 opens various channels of reactions, consisting of the neutron and 
a-captures. Because of much larger cross section, 17 soon captures another neutron to be 
14 C, and the latter is converted into 18 and then into 22 Ne by successive a-captures. The 
abundances of 17 0, 14 C, and 18 increase as 13 C is added, and yet, remain much smaller than 
the amount, AY" 13imix , of the mixed 13 C. In this case, 14 C attains to a fairly large abundance 
of ~ 0.2AYi3 jm i x since the rate of a-capture decreases relative to 3a-reaction rates with in- 
creasing temperature, as seen from Figure 4. At the same time, 17 reacts via the (a,n) 
channel to produce 20 Ne and releases neutron. This reaction is slower by a factor of ~ 10 4 
than the 13 C(a,n) 16 reaction, but, sufficiently faster than the 3a reactions. The released 
neutrons are mostly captured by 12 C and enter into the doubly neutron- recycling reactions of 
12 C(n,7) 13 C(a,n) 16 O(n,7) 17 O(a,n) 20 Ne. In addition to the neutron capture of 17 0, this cycle 
has a leakage of the 17 0(a,7) channel to 21 Ne, and it captures neutron to be 22 Ne with a fairly 
large cross section. 

The two paths from 17 meet at 22 Ne. The latter has a small cross section of neutron 
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capture, comparable with C and O, and the a-capture rates are slower than the 3a reactions 
for the temperature of T max ^3x 10 8 K. Accordingly, the flow of nucleosynthesis is bottlenecked 
and neutrons stagnate in 22 Ne. As a corollary, the abundance of 22 Ne reaches up to a half 
amount of the mixed 13 C and Y22 — AYi3. m i x /2 since two neutrons are spent to form 22 Ne. 

As neon isotopes accumulate, sodium and Mg isotopes are created, though in much 
smaller abundances than 22 Ne. There forms a steep drop in the abundances between 22 Ne and 
23 Na with the ratio V23/Y22 < <7 n7 ( 22 Ne)/cr n7 ( 23 Na) = 0.028, where the equality holds for a large 
neutron exposure of r ^> l/cr„ 7 ( 23 Na) = 0.48 mb -1 . Among Mg isotopes, the heaviest isotope 
Mg 26 tends to have larger abundance than 23 Na because of much smaller neutron capture cross 
section. The a-captures of Ne isotopes may produce the Mg isotopes directly by skipping 23 Na, 
though they work only partially at the temperature of this model. This also results in the 
largest abundance of 26 Mg because of larger rates and larger abundance of 22 Ne than those of 
20 Ne and 21 Ne. 

The elements heavier than Mg are formed solely by neutron captures (see Figure 4). They 
are produced late and the resultant abundances are much smaller than those of Na and Mg 
isotopes. Because of very small neutron capture cross section of 26 Mg, there forms another sharp 
drop of the abundances between 26 Mg and 27 A1 by more than Y27/Y26 < a n7 ( 26 Mg)/cr n7 ( 27 Al) = 
0.034. Among silicon isotope, 28 Si reaches as abundant as 27 Al and so does 31 P, reflecting their 
neutron capture cross sections much smaller than those of the ambient nuclei. Among the 
sulfur isotopes, 33 S is depressed because of the turn-around cycle via 33 S(n,a) 30 Si reaction. In 
comparison, 34 S has much larger abundance because of small neutron capture cross section. 

As soon as the mixing of 13 C ceases, the neutron density drops sharply to retard the 
neutron capture reactions, while the a-captures of the remainders of neutron-rich C and O 
isotopes keep adding the abundances of neon isotopes. In Models 4 and 5 of relatively high 
temperatures, neon isotopes also burn to form Mg isotopes. Accordingly, neutron supply may 
continue by 17 O(a,n) 20 Ne and 22 Ne(a,n) 25 Mg, although the neutron density is by far smaller 
and further decreases with the drop of temperature. 

Figure 6 shows the dependence of yields on the degree of mixing for two cases of Model 
1 (top panel) and Model 5 (bottom panel), with the lowest and highest temperatures, respec- 
tively. The hydrogen mixing makes striking difference in the variety and amount of synthesized 
elements when compared with the yields of the helium shell-flashes without the 13 C mixing, 
which is also plotted in this figure. Without the hydrogen mixing only the (0,7) reactions 
produce 16 0, and also, traces of 20 Ne and 24 Mg for the flashes of high temperature. In contrast, 
the hydrogen mixing enables to produce not only the heavier a nuclei in much larger amounts, 
but also, their neutron-rich isotopes and the odd nuclei appreciably. 

The amounts of yields increase with the amount of mixed 13 C. There is a general ten- 
dency that the dependences grow stronger for heavier elements, though varying with elements. 
Among the oxygen isotopes, 16 O is produced via the recycling reactions, and the resultant 
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abundance depends rather weakly on the mixed amount of C. There is little difference among 
the models. On the other hand, the heavier isotopes differ in their abundances greatly between 
Model 1 and Model 5 since they are mostly destroyed by a-captures after the cease of mixing. 
The latter is also the case for nitrogen isotopes. In contrast, the fluorine depends weakly on 
the models and increases largely with the mixing amount of 13 C, which may result from the 
competition of two formation paths (see below). 

Among the neon isotopes, the heaviest isotope 22 Ne increases nearly in proportion to 
the amount of mixed 13 C, as stated above. In comparison, the lighter isotopes, 20 Ne and 21 Ne, 
have weaker dependence on the amount of mixed 13 C as a result of the competition between 
the a and neutron captures of 17 0; for a smaller mixing of 13 C, the a-capture dominates over 
the neutron capture to produce more 20 Ne than 22 Ne and vice versa for a larger mixing. The 
temperature dependences are discernible in the larger abundances of 20 Ne for Model 5 and in 
the survival of neutron-rich carbon and oxygen isotopes and of nitrogen in larger amounts for 
Model 1. 

Beyond neon, the dependence on the amount of mixed 13 C grows steeper with the mass 
number while the relative distribution abundances approach to the steady-state distributions of 
neutron captures for larger mixing. Among the magnesium isotopes, the heaviest isotope 26 Mg 
has the largest abundance, and the influence of a-captures is discernible for the shell flashes of 
the higher peak temperatures. For still heavier elements, the dependence remains very strong 
even for the largest mixing. The abundances of phosphorus and sulfur isotopes vary by more 
than a third power of the amount of mixed 13 C. 

Figure 7 compares the yield abundances, resulting from different mixing durations for 
the same amount of mixed 13 C. The abundance patterns are almost independent of mixing 
duration for sufficiently long mixing epoch, while the production of heavier nuclei is suppressed 
for the shorter durations of At ^ 10 7 s for Model 5 and 10 10 s for Model 1. These durations are 
related to the timescales necessary for a-captures to synthesize 22 Ne that acts as seed for the 
synthesis of heavier nuclei. This fact manifests itself in the smaller abundance of 20 Ne for the 
short mixing durations. Further, the mixing duration also has an effect on the a -captures after 
the end of 13 C mixing. In Model 5, 25 Mg and 26 Mg are more abundant for shorter durations 
since they are formed via a-captures of 22 Ne after the neutron irradiation has ceased. In Model 
1, we see that 14 C or 14 N, 15 N, 17 and 18 remain with larger abundances for longer durations, 
which is simply because of shortage of time and temperature to burn them via a-capture after 
the cease of mixing. 

Figure 8 shows the variations of neutron exposure during the shell flashes. It is deter- 
mined essentially by the amount of mixed 13 C and increases nearly in proportion to its square 
root because of concomitant production of neutron absorbers. The weak dependence on the 
mixing duration results from the competition between the a and neutron captures of 17 0. For 
the shorter durations, large neutron density suppresses the second-recycling reaction while for 
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the longer duration, the production of heavier elements increases the average cross sections of 
neutron capture, both of which work to diminish the neutron exposure. 

In summary, in the nucleosynthesis during the helium flashes with the hydrogen mixing, 
neutrons play a critical role in accelerating a-captures to create the a nuclei of oxygen through 
magnesium as catalytic agents, and also, to generate their neutron-rich isotopes and odd nuclei 
as incorporated. In the following, we discuss the formation paths to the light a-nuclei, their 
neutron-rich isotopes and related odd nuclei separately. 

3.1. Formation Path to Particular Elements 

3.1.1. Oxygen Isotopes and Fluorine 

To create 16 0, there are two paths, (0,7) and (71,7) plus (a,n), from 12 C in the helium 
zone with the hydrogen mixing. The 13 C(a,n) 16 reaction is faster by many orders of magnitude 
than the 12 C(a,7) 16 (see Figure 4), and the ratio of production rates is estimated at A? 6o — 10 5 
for the Model 4 of A( 13 C/ 12 C) = 1CT 2 in Fig 5. Accordingly, the (71,7) and (a,n) reactions are 
the main channel to 16 at least while the 13 C mixing continues. 

The amount of 16 O, thus produced, results from the competition between the neutron 
capture by 12 C and 16 O. Neutrons, if captured by 12 C, enter into the neutron-recycling reaction 
12 C(n,7) 13 C(a,n) 16 0. On the other hand, if absorbed by 16 0, they escape from the recycling 
circuit and finally end by settling in the neutron-rich isotopes from 14 C up to 22 Ne. By neglecting 
a small variation of 12 C abundance, we may evaluate this competition analytically (see Appendix 
1). For r < l/cr n7 ( 16 0) = 26 mb" 1 , the resultant oxygen abundance is given as a function of 
the amount of mixed 13 C from eq. (A14) by 

1/2 

(Ay 13)mix /^i2) 1/2 , (9) 

where £ denotes the ratio of reaction rates between the a to neutron captures of 17 and / 
denotes the branching ratio of the (0,7) to the total reactions with a particles. For still larger 
r, an upper bound is set at a steady state value of 

Y m /Y 12 < a n7 ( 12 C)/a„ 7 ( 16 0) = 0.40. (10) 

On the other hand, the 12 C(a,7) 16 reaction works in longer time, producing an amount 
of SY 16 = 2.9 x 10~ 3 in total and <5Y"i 6 = 2.19 x 10~ 3 after the mixing has ceased for Model 4 in 
Figure 5. As the convection extends outwards during the helium flashes, carbon and oxygen 
suffer from dilution with helium in the helium zone newly incorporated. Consequently, as the 
shell flashes recur without hydrogen mixing, the ratio of carbon to oxygen approaches their 
production rates, C and O, due to the 3a and 12 C(a,7) 16 reactions, respectively (Suda et al., 
2004). Thus, we have 

Y 16 /Y 12 ~ 6/C ~ (0.07-0.4)0.25(A 12 /0.2)(0.8/F) 2 (10 3 g cm" 3 //)) (11) 

for the relevant temperature range (logT = 8.2 ~ 8.5) where Y and A 12 are the helium and 
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2(1+0 ^ 7 ( 12 C) 
2 + £(l + /)a n7 (i60) 



carbon abundances in the helium convective zone [Y + X\2(l + O/0.75C) = 1]. This sets the 
range of oxygen abundance resulting from the third dredge-up only. 

Other isotopes, 17 and 18 0, are produced by the 16 0(n,7) 17 reaction and by the 
subsequent reactions of 17 0(n,a) 14 C(a,7) 18 0, respectively. Because of much larger cross section 
of neutron capture, 17 remains much smaller in the abundance than 16 [Y17/Y16 < 0.0038/(1 + 
£)] while irradiated by neutrons. After the cease of 13 C mixing, it is effectively destroyed by the 
a-captures. The same is true for 14 C (or its daughter nuclei 14 N) and 18 0. They accumulate 
as 13 C is added, and decreases rapidly after the neutron supply stops. The abundance of 18 
finally reduces much smaller than that of 17 0, which reflects much larger a-capture rates of 
18 0. 

As 18 increases, 19 F is formed, though much smaller in the abundance. Under neutron- 
rich circumstance, there are two channels of 18 0(n,7) 19 0(e~z/) 19 F and 15 N(a,7) 19 F (see Mowlavi 
et al., 1998). The latter a-capture channel dominates the former channel even for large 13 C 
mixing; the latter contribution accounts for 2/3 ~ 0.9. This can be discernible in Figure 5, 
where the abundance is a few times larger than expected from that of 18 O under the steady- 
state of neutron capture during the mixing and increases rapidly via a-capture of 15 N after the 
cease of mixing and the resultant drop of neutron density. Although the abundance is itself not 
large, this can give rise to large enrichments because of its smaller solar abundance, reaching 
the enrichment comparable with carbon for the largest mixing. 
3.1.2. Neon isotopes and Sodium 

The argument similar to the 16 case is applicable to the effective paths to create 20 Ne. 
Since A2o Ne ~ 10 4 in Model 4 of Fig 5, the main path is expected to be the (n,7) and (a,n) 
reactions as well. Differently from 16 case, however, the a-capture of 17 has the (0,7) 
branch, which works as effective neutron leakage out of doubly neutron- recycling circuit. Thus, 
the abundance ratio is determined by the branching ratio as Y21/I20 — 1/5-5 — 1/10 around 
the temperatures (T§ ~ 2 — 3) of interest here (Caughlan & Fowler, 1988; Angulo et al., 1999). 
In addition, the formed 21 Ne may absorb neutron while 20 Ne largely survives for the large 
neutron exposure. As neutron capture proceeds, the abundance ratio decreases and approaches 
the steady-state values of Y^i/Y^q — 0.08/(1 — /), which happens to be only slightly smaller 
than given by the branching ratio. It also decreases slightly at high temperatures because of 
larger 21 Ne(a,n) rate than 20 Ne(a,7) rate. The abundance ratio between 20 Ne and 22 Ne results 
essentially from the competition between the (a,n) and (n, a) reactions of 17 0, and hence, 
Y20/Y22 = (1 - /)£■ Then 22 Ne overweighs 20 Ne except for small mixing of A( 13 C/ 12 C) < 0.001 
since £ decreases for greater mixing, and attains at the abundance of 

Y 22 ^ AY 13Mx /2, or [Ne/C] ~ 0.21 + log A( 13 C/ 12 C) (12) 

unless the peak temperature reaches high enough to promote a-captures of 22 Ne. Here we adopt 
the solar abundance of Asplund et al. (2005). 
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Since Na is formed by neutron capture, the abundance is limited by the ratio of neutron 
cross sections, and from eq. (12) we have 

Y 23 < 0.0 281- 22 ~ 0.014AF 13imix , or [Na/C] < 0.33 + logA( 13 C/ 12 C), (13) 

which may be compared with the enhancement during the third dredge-up in eq. (A17) of 
Appendix 2. 

3.1.3. Magnesium isotopes and Aluminum 

As for the creation of 24 Mg, the difference between the 20 Ne(a,7) and 21 Ne(a,n) reactions 
is not so large, giving A24 Mg — 30 — 700 in Fig 5. Furthermore, these reactions are slower than the 
3a reaction for logT < 8.45 (see Figure 4). For weak flashes with such low peak temperatures 
as in Model 1, therefore, the dominant process to create 24 Mg is the neutron capture of 23 Na 
and the subsequent /3-decay. The magnesium isotopes keep increasing with Na, and eventually, 
26 Mg becomes more abundant than 23 Na because of much smaller neutron capture cross section. 
Under the steady-state condition, the ratio of magnesium to sodium may reach 

Y 26 /Y 23 < a n7 ( 23 Na)/a n7 ( 26 Mg) ~ 16.7, or [Mg/Na] < -0.12. (14) 

Here [Mg/Na] includes the contributions of all stable Mg isotopes. For strong shell flashes 
with such high peak temperatures as in Model 5, the a-captures may compete or overweigh 
the neutron captures, in particular when the mixing is weak and the neutron density is small. 
Because of larger a-capture rates of 22 Ne than that of the other neon isotopes, 25 Mg and 26 Mg 
are produced more than 24 Mg. In this case, the 26 Mg has also the largest abundance among 
the magnesium isotopes, since 25 Mg captures neutron with much larger cross section. But 
note that a-capture reactions of 22 Ne are subject to large uncertainties about three orders of 
magnitude in the largest, as seen from the NACRE compilation. Recently, the uncertainty has 
been reduced by the reanalysis but still remains as large as a factor of ten (Koehler, 2002). 

Consequently, the peculiar abundance ratios among the magnesium isotopes result from 
the difference in neutron capture cross sections as well as from the contribution of a-capture 
reactions of 22 Ne. The most abundant isotope is 26 Mg, while the other isotopes, 24 Mg and 
25 Mg, turn out to be as abundant as 23 Na, reflecting the nearly similar neutron capture cross 
sections. Because of the very small cross section of 26 Mg, an additional steep decrease in the 
abundances results between 26 Mg and 27 Al. The upper limit of the ratio is given by 

Y 27 /Y 26 < a„ 7 ( 26 Mg)/cr n7 ( 27 Al) ~ 1/30, or [Al/Mg] < -0.34, (15) 

where the contributions from all Mg isotopes are included in the [Al/Mg] ratio. 
3.1.4- Si isotopes and heavier elements 

Beyond magnesium isotopes, the elements are created solely through a series of neutron 
captures and /3-decays since the a-capture reactions can no longer make any significant contri- 
bution under the temperature realized during helium shell flashes, as seen from Figure 4. The 
fact expresses itself in the abundance distributions of Al through S isotopes in Figures. 6 and 7. 

14 



All the models show a general trend that their abundances decrease for heavier nuclide and also 
for larger neutron capture cross sections. The abundances increase with the amount of mixed 
13 C more rapidly than linearly and approach the steady-state ones inversely proportional to 
the cross sections. Accordingly, the silicon isotopes are less abundant than, or comparable to, 
Al, and local peaks appear at phosphorus and 34 S, which have much smaller neutron capture 
cross sections as compared with ambient nuclides. 

4. Comparisons with the Observations of HMP/UMP and EMP Stars 

We discuss the relevance to the abundance patterns observed for the stars in the Galactic 
halo from the viewpoint of the binary scenario (Suda et al., 2004). In this scenario, the nuclear 
products in the helium zone are eventually carried out to the surface of erstwhile primary 
st result of the deep mixing and dredge- up (He-FDDM), which occurs after the shell 

flash grows strong enough to engulf hydrogen at sufficiently high rate. In addition, the nuclear 
products are also dredged up to the surface during the helium shell-flashes (TDU) subsequent 
to He-FDDM in some cases. Then, the matter in the envelope of the primary stars is planted 
on the low-mass secondary members through wind accretion to be observed as the surface 
pollution. The yields may suffer dilution with the envelope matter both in the convective 
envelope of primary stars when they are dredged up and in the surface convection of secondary 
stars when they are accreted. Accordingly, it makes sense to compare the relative abundances 
among the elements. 

Figure 9 shows the abundance patterns of yields of our computations for various amounts 
of mixing between A( 13 C/ 12 C) = 10~ 4 and 0.1 (solid lines). Note that the elemental abundances 
are plotted so as to give the same carbon enrichment relative to iron as observed for HE1327- 
2326 ([C/Fe] =4.26). In addition to carbon, oxygen is also enhanced as large as [0/C]~— 2~— 1 
with rather weak dependence on the amount of mixed 13 C. The neon enhancement is in 
proportion to the amount of mixed 13 C. For heavier elements, the enhancement decreases in 
general and the dependence on the mixing amount grows stronger. Sodium, magnesium, and 
aluminum are all enhanced for A( 13 C/ 12 C) > 0.01 while only sodium for A( 13 C/ 12 C) ~ 0.001. 
Fluorine and phosphorus are greatly enriched for A( 13 C/ 12 C) > 0.01, reaching the comparable 
enhancement as carbon for A( 13 C/ 12 C) ~ 0.1 because of their small solar abundances. 

In comparing the yields from helium shell flashes with the observed abundances, we 
should take into account some later modifications that may be introduced when and after 
they are dredged up to the surface. During He-FDDM, a hydrogen shell flash ignites in the 
middle of helium convective zone, and the largest effect is to burn carbon into nitrogen as much 
as [N/C] ~ 0.7 ~ —0.3 in the hydrogen-flash convection (Fujimoto et al., 2000; Iwamoto et al., 
2004), while other heavier elements are mostly left unaltered because of larger Coulomb barriers 
and a short duration of hydrogen flash-convection. Hydrogen mixing and He-FDDM works for 
stars of mass M ^ 3M Q until the CNO abundances in the envelope exceeds [CNO/Fe] ^ —2.5 
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(Fujimoto et al, 2000; Suda & Fujimoto, 2008). For stars of mass M > 1.5 M , TDU follows 
He-FDDM to enrich the envelope with 12 C and 16 but without nitrogen. In the envelope of 
massive AGB stars, the nuclear products, after dredged up to the surface, are further processed 
by hydrogen burning in the envelope convection (hot bottom burning; HBB). There are the 
lower limits of stellar mass for TDU and for HBB, and are argued to decrease for smaller 
metallicity. It is true, however, that they are subject to uncertainties due to the treatment of 
convective mixing, such as the mixing length theory and overshooting. We rather consider that 
these parameters should be determined from the comparison with the observations of EMP 
stars. 

The observed abundances for three HMP/UMP stars are also plotted by symbols. For 
HE1327-2326, the pattern of abundance enhancement from carbon through aluminum is well 
reproduced by our yields with the amount of mixing A( 13 C/ 12 C) ~ 0.02 as long as the LTE 
abundances are concerned (filled circles). Nitrogen enrichment ([N/C] =0.3; Aoki et al. 2006)is 
also consistent with the processing by hydrogen shell-flash during He-FDDM, as shown by 
Fujimoto et al. (2000) for 0.8 M Q with Z = and [Fe/H] = -4, and by Iwamoto et al. (2004) 
for 1 — 2M and [Fe/H] = —2.7, if we take into account the differences in the metallicity and 
the core mass. The non-LTE correction (open circles) gives smaller abundance to sodium while 
larger abundance to aluminum than the LTE abundances both by ~ 0.4 dex, but remains within 
the 2a of typical errors. 

For HE0107-5240, the operation of TDU subsequent upon He-FDDM is proposed by 
Suda et al. (2004) to explain the small nitrogen abundance by a factor of ~ 30. As the surface 
carbon abundance increases ~ 30 times by TDU, the surface nitrogen abundance is reduced as 
much but the surface abundances of the other yields are reduced by a factor of ~ 20 (dashed 
line). The difference in the reduction factors is related to the fact that without hydrogen 
mixing, nitrogen is burnt out but the other yields can survive in the convective zone during the 
subsequent helium shell-flashes to be dredged up at the same time with carbon (see Suda et al., 
2004). The effect of dilution by TDU is consistent with the observed sodium enhancement. 
The magnesium abundance with rather small enhancement and an upper limit of aluminum 
abundance are also consistent with the dilution. 

As for HE0557-4840, the enhancement is observed only for carbon ([C/Fe] = 1.7 — 2.1), 
but the excess is much smaller as compared with two HMP stars; nitrogen and oxygen are given 
only the upper limits. There are no indications for the enrichment for sodium, magnesium or 
aluminum (Norris et al., 2007). The observed abundance pattern does not necessarily demand 
the hydrogen mixing events and is explicable in terms of TDU alone. 

In the following subsections, we discuss the comparisons on the carbon and oxygen ratios 
and the abundances of light elements of sodium, magnesium and aluminum in some more detail 
for the HMP/UMP stars, and also for EMP stars. In our computations, we assume the metal- 
free stars, and yet, the results may be little affected by the presence of pristine metals as long 
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as their contributions to neutron captures are much smaller than that of 0. Since the pristine 
CNO elements are converted into 22 Ne in the helium convective zone, the latter condition is 
reduce to 

Y 56 < Vi 6 ,o^n 7 ( 16 0)/a n7 ( 56 Fe) ~ O.O27F 5 6,0(n6,o/O.Oir 12 )(X 12iHc -conv/O.2). 

Accordingly, the present results may be applicable to EMP stars of [Fe/H] < —2.5, and we take 
them into account. On the basis, we finally discuss the origin of three HMP/UMP stars. . 

4-1. carbon- oxygen enrichment 

Our results suggest that the helium shell flash with hydrogen mixing is an effective site 
for oxygen production. The neutron-recycling reactions can bring the oxygen and carbon ratio 
up to eq. (10), and hence, 

[O/C] < -0.66, (16) 

with the use of the solar abundances by Asplund et al. (2005), although it demands larger 
neutron exposure, or larger mixing of 13 C, than in our computation to realize the ratio close 
to it. This upper bound is larger by a factor of ~ 10 than the ratio in eq. (11) that can be 
attained by TDU alone; 

[O/C] = (-2.5-- 1.75) +log[(X 12iHc -conv/0.2)(0.8/F H c-conv) 2 /(p/3 x 10 3 g cm" 3 )], (17) 

where -X"i2,He-conv and Yke-conv denote the abundances of carbon and helium in the helium flash 
convective zone, respectively, and the density is taken from the middle stages of the strongest 
shell flashes (Models 4 and 5) since TDU occurs for relatively massive stars. 

Figure 10 shows the comparison of these theoretical predictions from the AGB nucle- 
osynthesis and dredge-up with the observed ratios of oxygen to carbon for HMP/UMP and 
EMP stars. The upper bound for He-FDDM in eq. (16) and the upper and lower limits for 
TDU in eq. (17) are drawn by thick solid and broken lines, respectively. For HE1327-2326, 
we take the carbon abundance of [C/Fe] = 4.26 ±0.24 from the ID LTE analysis of CH lines 
(Aoki et al., 2006). For the oxygen abundances, an upper bound of [O/Fe] < 3.1 is derived 
from non-detection of both OI triplet at A7771 — 5 and [OI] at A6300, while the abundances 
[O/Fe] = 3.7±0.2 and 2.8 ±0.2 result from the ID LTE analysis of UV OH features and from 
the applications of 3D corrections, respectively (Frebel et al., 2006). These ratios are consistent 
with the range that can be reproduced by the hydrogen mixing and dredge-up by He-FDDM. 
Even the largest one is consistent with our computation within an error if we regard nitrogen 
as carbon processed by hydrogen shell-flash during He-FDDM . The ratios, including the upper 
limit for OI line, lies well above the range for TDU. For HE0107-5240, we have the carbon 
and oxygen abundances only from the ID LTE analyses at [C/Fe] = 4.0 ±0.2 (Christlieb et al., 
2004), and [O/Fe] = 2.3 ± 0.2 from UV OH lines (Bessell et al, 2004), which gives slightly 
smaller ratio of [O/C] = —1.7 ±0.3 than for HE1327-2326. This ratio lies just below the upper 
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border of the range in eq. (17), and is explicable in terms of TDU, if allowance is made for 
errors in the abundance determination. 

For EMP stars, we plot all the samples that we find the measured oxygen abundances 
in the literature with use of SAGA Database (The Stellar Abundance for Galactic Archeology 
Database; Suda et al., 2008). There is a general trend that the O/C ratios decrease with carbon 
enrichment. The largest ratios of [O/C] ~ 0.5 — 1.5 are found for EMP stars without the carbon 
enhancement ([C/Fe] < 0.5), among which the oxygen abundances are in the range between 
[O/Fe] ~ 0.5 and 1.0 and the variations in the [O/C] ratios are attributable to the conversion 
of carbon into nitrogen possibly by CN cycle processing (Spite et al., 2005). 

A salient feature of the oxygen abundances in the EMP stars is that such large enrich- 
ments as exceeding [O/Fe] ~ 1 are observed only for the carbon-enhanced EMP (CEMP) stars 
of [C/Fe] > 0.5. Among the EMP giants without carbon enhancement, the largest abundance 
is [O/Fe] = 1.14 with the largest ratio [O/C] = 1.41 for BPS BS 16934-002 in this figure (Aoki 
et al., 2007, see also SAGA Database). The CEMP stars exhibit lower O/C ratios, and most 
of them have such low ratios as explicable in terms of the AGB nucleosynthesis with hydrogen 
mixing and He-FDDM and/or with TDU. In this figure, we delineate the ranges of [O/C] for 
the various mixtures of the matter, dredged up by He-FDDM or by TDU with the envelop mat- 
ter of the pristine oxygen abundance taken to be [O/Fe] = 0.5 or 1.0, as a function of carbon 
enrichment. CEMP stars divide into two subgroups according to the s-process enrichment, i.e., 
CEMP-s with [Ba/Fe] > 1.0 (filled symbols) and CEMP-nos with [Ba/Fe] < 1.0 (open sym- 
bols). All the CEMP-s stars fall in the domain of O/C ratios realized by He-FDDM and TDU, 
mostly in the range that cant be attained only by He-FDDM, and about a half of CEMP-nos 
stars also exhibit small O/C abundances consistent with the AGB nucleosynthesis. Two of the 
smallest O/C ratios (HE0024-2523, CS22958-042) are quite low values reproducible by TDU 
alone. HE0024-2523 belongs to a binary of period 3.14 days, indicative that it has under- 
gone the Roche-lobe overflow and the common envelope phase (Lucatello et al., 2003), while 
CS22958-042 is a CEMP-nos star but shows the largest sodium enrichment (Sivarani et al., 
2006). In addition, for CS30322-023, which has a large nitrogen enrichment relative to carbon 
([N/C] = 2.21; Masseron et al. 2006), the O/C ratio reduces as small as consistent with the 
TDU nucleosynthesis without hydrogen mixing if nitrogen is regarded as having been produced 
from carbon ([0/C(+N)] ~ —1.7). If this is the case, the observed large ratio of nitrogen to 
carbon demands a deep processing by CN cycle, and the most likely site is the hot bottom 
burning in the envelope of massive AGB stars. This in turn supports the binary scenario along 
with the variations in the radial velocity of period 191.7 days, at variance with the argument 
by Masseron et al. (2006) that this star is a currently TP- AGB stars of mass ~ 0.8 M Q . 

On the other hand, some CEMP-nos stars show the O/C ratios and the enrichment of 
O/Fe too large to be explicable in terms of the AGB nucleosynthesis, and seem to comprise 
the groups that demand other oxygen source(s). Among them are two giants, CS29498-043 
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and CS22949-037, both greatly enriched with nitrogen as [N/Fe] = 2.27 (CS29498-043; Aoki 
et al., 2004) and 2.59 (CS22949-037; Depagne et al, 2002) as well as with Na and Mg greatly 
but not with Al (see below). If nitrogen is added to carbon, the O/C ratio prior to the 
processing by CN cycles may decrease and yet, are still too large to be reproduced by the AGB 
nucleosynthesis. Another CEMP-nos subgiant, HE1300+0157, exhibits a much larger oxygen 
enrichment, [O/C] = 0.4 ([O/Fe] = 1.8) than expected from the AGB nucleosynthesis, although 
the overabundances may be alleviated if 3D effect is included in the abundance analysis (Frebel 
et al., 2007). This source is different from the above two in the absence of the enrichment of Na- 
Al. In addition, a CEMP-nos dwarf, CS 31080-095, also exhibits a slightly larger enrichment 
of oxygen relative to carbon than expected from our results. It may be distinct from the other 
three sources in larger metallicity by ~ 1 dex and also in a moderate enrichment of neutron 
capture elements ([Ba/Fe] = 0.77). 

In summary, the AGB nucleosynthesis with the hydrogen mixing in combination of 
pristine oxygen abundance of [O/Fe] ~ — 1 can provide reasonable explanation to the observed 
carbon and oxygen abundances for the majority of CEMP stars (all the CEMP-s stars and half 
of CEMP-nos stars) including HMP stars. In addition, our results may serve as a filter to single 
out such stars as demand other sources of carbon and oxygen, and identify the nucleosynthetic 
signatures presumably of Pop. Ill and early-generation supernovae. 

4-2. Na, Mg and Al enrichment 

Among the light elements that can be enriched by the neutron-recycling reactions, a 
number of observations exist only for sodium, magnesium and aluminum. Figure 11 shows a 
summary of the abundances of these elements, obtained by our numerical computations, and the 
comparison with the observations of HMP/UMP and EMP stars in the [Mg/Na] and [Al/Mg] 
diagram. The numerical results (open circles) fall in the left-bottom quarter demarcated by 
vertical and horizontal lines, which denote the steady-state ratios under the neutron captures, 
respectively ([Mg/Na] = -0.12 and [Al/Mg] = -0.34). 

In this figure, we also delineate the evolutionary paths of nucleosynthesis for several 
models by solid curves. Those curves which run most leftwards represent the evolutionary paths 
for which neutron captures dominate over a-captures in producing magnesium isotopes, as seen 
for Model 1 with A( 13 C/ 12 C) = 0.1 and At mix = 10 12 s and for Model 5 with A( 13 C/ 12 C) = 0.1 
and At mix = 10 6 s. Along them, the ratios of Mg/Na and Al/Mg both increase and approach 
the steady-state values at the intersection of vertical and horizontal lines. The approach to the 
steady-state abundance ratio is faster between Al and Mg than between Mg and Na since the 
ratio of the neutron capture cross sections of Al to Mg is much larger than that of Mg to Na, 
and hence, much smaller abundance ratio is sufficient for the former than for the latter to reach 
the steady state. The paths delineate the left boundary of the distribution of our yields which 
we call "the neutron capture path" in the following. 
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In comparison, some models run on the right side of the neutron-capture path, as seen 
for Model 5 with A( 13 C/ 12 C) = 0.1 and At mix = 10 9 s and for Model 4 with A( 13 C/ 12 C) = 0.1 
and At m ; x = 10 10 s. This is the case for the models with high peak temperatures and long 
mixing duration (or small neutron density). Higher temperature and smaller neutron density 
promote the 17 0(a,7) 21 Ne(a,n) 24 Mg reactions to make the contribution to the production of 
24 Mg larger. Thus the models trace courses with larger [Mg/Na] for given [Al/Mg] in this 
figure. As 22 Ne accumulate and the neutron captures with 22 Ne as seed come to dominate over 
the Mg production, our models approach to the neutron-capture path and eventually join in it. 

After the cease of mixing of 13 C, the a-capture reactions proceed if the temperature 
remains sufficiently high. For the models of short durations of mixing with the high peak 
temperatures, therefore, the evolutionary path departs again from the neutron-capture path 
and runs to right- and down-ward as magnesium isotopes are formed by a-captures of 22 Ne. 
The shorter the duration of mixing, the greater is the departure of the end products from the 
neutron-capture path. Two dash-dotted lines denote the resultant abundance ratios for Model 
5 (upper line) and Model 3 (lower line) both with A( 13 C/ 12 C) = 0.01. 

The HMP stars exhibit the abundances of Na, Mg and Al well within the range that can 
be reproduced by the neutron-capture reactions, which are quite different from the averaged 
values that are observed for EMP stars without the carbon enrichment ([Mg/Na] ~ +0.5 and 
[Al/Mg] ~ — 1). HE1327-2326 is located right on the neutron capture path when we adopt 
LTE abundance ratios. HE0107-5240 falls near but right to HE1327-2326 despite that the 
enrichments of Na and Mg are smaller by ~ 1.6 dex and more. Since the Al abundance is 
an upper limit, the abundance ratio lies below the neutron-capture path, which is taken to 
be indicative of the shell flash of higher temperatures. On the other hand, the UMP star, 
HE0557-4840, displays the absence of significant enrichments for any of these elements and is 
located in the middle of the EMP stars without carbon enrichment. 

The distribution of EMP stars on this diagram is clearly divided between CEMP-s 
and CEMP-nos stars. All the CEMP-s stars fall within, or near to, the ranges of Mg/Na 
and Al/Mg ratios, encompassed by the numerical results, except for the two stars of short- 
period binaries, HE0024-2523 (period of 3.14days Lucatello et al., 2003) and CS30322-023 
(possibly 192days Masseron et al., 2006), which are thought to have undergone the common- 
envelope phase. This is consistent with the fact that the s-process nucleosynthesis is a necessary 
consequence of neutron capture reactions under the existence of pristine seed nuclei. The 
observed values of [Mg/Na] are shifted slightly rightwards of larger magnesium abundance 
from the neutron capture path. This may be attributable to the contribution of pristine 24 Mg 
which has presumably larger abundance than sodium. In addition, most of CEMP-s stars 
have the sodium and carbon ratios similar to, or below, the value [Na/C] ~ —1.7 observed 
for HE1327-2326. Therefore, the origin of these elements is explicable in terms of hydrogen 
mixing of A( 13 C/ 12 C) < 0.02, while two stars, LP625-44 and CS22942-019, have larger sodium 
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to carbon ratios ([Na/C] = -0.50 and -0.78). 

In contrast, CEMP-nos stars distribute on the bottom-right quarter in this figure with 
smaller Na abundance relative to Mg than predicted from the neutron- capture nucleosynthesis; 
most have the abundance ratios of [Al/Mg] and [Mg/Na] similar to EMP stars without the 
carbon enhancement, but with larger scatters. The two CEMP-nos stars with large O/C ratios, 
CS29498-043 and CS22949-037, discussed above, are characterized by the large enhancement 
of Na and Mg. Depagne et al. (2002) and Aoki et al. (2004) discuss the peculiar abundance 
patterns of these stars in connection with a peculiar supernovae such as proposed by Umeda 
& Nomoto (2003) and also with the rotational effects as asserted by Meynet & Maeder (2002). 
These stars have the Al/Na ratios similar to the CEMP-s stars but the magnesium abundances 
are larger by ~ 1 dex. The comparison with the evolutionary paths of Model 5 in the figure 
suggests the nucleosynthesis during helium shell flashes of larger core masses with smaller 
neutron density for which the temperatures reach high enough to promote the a-capture of 
22 Ne. Komiya et al. (2007) propose that CEMP-nos stars are formed from the binaries with 
a primary star of mass M > 3 M and also with still more massive primary that develops an 
ONe core and ends with a supernova explosion. The binary origin with the massive primary 
can be an alternative, while we have to seek for the source of large oxygen enrichment relative 
to the carbon enrichment. The abundance patterns of these stars wait for further study, both 
theoretical and observational, which may shed light on the new aspect of the star formation 
and chemical evolution in the dawn of the universe. 

It is worth noting that the stars, CS22958-042, CS30322-023, and HE0024-2523, with 
such low 0/C(+N) ratios as consistent with TDU alone, display quite different behaviors in 
this diagram from each other; the largest enrichment of sodium ever known for CS22958-042 
(Sivarani et al., 2006), large enrichment both of aluminum and sodium for CS30322-023 and 
no enhancement either of sodium nor aluminum for HE0024-2523. While they share the ni- 
trogen enrichment ([N/Fe] = 2.15, 2.81, and 2.1), the nitrogen to carbon ratios differs greatly 
([N/C] = —1.0, 2.21 and —0.5). The enhancement of s-process elements is absent for CS22958- 
042, while very large for HE0024-2523 and CS30322-023 with large [Pb/Ba]( = 1.84,0.95), respec- 
tively. These may reflect the range of AGB nucleosynthesis according to the common-envelope 
evolution since these all are suspected as the members of short period binary systems from the 
radial velocity variations. 

Although we have discussed so far about the LTE abundances, it is pointed out that 
departures from LTE have a significant impact on the abundance determination by use of 
the resonant lines. For HE1327-2326, the NLTE corrections give the Al/Mg ratio larger 
than that of LTE abundances by ~ 0.4 dex and exceeding the upper limit reached via neu- 
tron capture reactions. In order to explain the NLTE ratio within the present framework of 
neutron capture reactions, we need as much smaller ratio of neutron capture cross sections, 
cr ri7 ( 27 Al)/o" n7 ( 26 Mg), which possibility is not excluded from the current uncertainties of their 
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cross sections. Otherwise, the enrichment of Al may be realized by proton captures with the 
burning of only a small fraction of Mg, but we have to seek for the appropriate site. It is 
thought, however, that the abundance determination is subject to systematic errors also due 
to the 3D effects and this may have opposite effects to NLTE corrections (Collet et al., 2006). 
We had better wait for future studies to consider the effects of departures from LTE and the 
uncertainties of nuclear reaction rates. 

4-3. Origin of Abundance Patterns of UMP/HMP Stars 

We may now discuss the constraints on the origin of HMP/UMP stars within the frame- 
work of the binary scenario (Suda et al., 2004), based both on the above discussion on the 
abundance patterns and on the general picture of evolution of metal-free and extremely metal- 
poor stars, presented by Fujimoto et al. (2000, see also Komiya et al. 2007 and Suda & 
Fujimoto 2008). 

For HE1327-2326, the observed abundance pattern of the light elements from carbon 
through aluminum can be well reproduced by the neutron-recycling reactions in the metal-free 
or extremely metal-poor AGB stars. The abundance ratios among Na, Mg and Al reflect the 
yields of neutron-capture reactions for the model of the lowest temperature with A( 13 C/ 12 C) ~ 
0.02, and nitrogen is produced during He-FDDM. The low-mass nature of the primary star is 
in accordance with the absence of the dilution due to TDU for these elements, which suggests 
that HE1327-2326 was formed as, and evolved from, the second member of a low-mass binary 
system with a primary star of mass 1.0 M Q < M < 1.5 M Q . For such a low- mass primary, it is 
possible to enrich the envelope with carbon up to, or even beyond, the observed enrichment 
of [C/H] = -1.4 by He-FDDM without TDU (Fujimoto et al, 2000; Iwamoto et al, 2004). In 
addition, our interpretation is consistent with the observed depletion of Li since Li is readily 
destroyed during the AGB phase of primary star even if it has been formed during a hydrogen 
shell-flash of He-FDDM; The s-process nucleosynthesis and Sr enrichment will be discussed in 
a subsequent paper. 

For HE0107-5240 abundances (Christlieb et al., 2004), the large enrichment of carbon 
([C/H] = —1.28) and the sodium enrichment ([Na/Fe] = 0.81) along with the small nitrogen 
enrichment [N/C] = -1.72 to -1.13 may require both He-FDDM and TDU (Suda et al., 2004). 
Accordingly, HE0107-5240 evolved from a binary system with the primary of mass between 
1.5 M Q < M < 3M . The oxygen abundance is explicable either by the hydrogen mixing or by 
TDU. The sodium enrichment itself is relatively small in relation to carbon ([Na/C] = —3.19) 
and explicable by TDU alone as seen from eq. (13). With such a small pristine abundance 
([Fe/H] = —5.3 with normal CNO abundance), however, TDU cannot grow strong enough to 
enhance the carbon abundance up to the observed value (Suda & Fujimoto, 2008). Accordingly, 
it must have experienced the hydrogen mixing and resultant He-FDDM to increase the carbon 
abundance above [C/H] > —2.5. Suda et al. (2004) predict the binary period of this system 
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at P or b = 76-150 yrs, and for HE1327-2326, the binary period can be even longer because of 
shallower surface convective zone (Komiya et al., 2007). 

As for HE0557-4840, few information is available except for the relatively small carbon- 
enrichment of [C/H] = —3.1 — 2.63, smaller by an order of magnitude than the other two stars. 
The upper limit of nitrogen abundance excludes the low- mass primary star with He-FDDM since 
the ratio of N/C has to be significantly smaller than observed from HE1327-2326. Differently 
from HE0107-5240, the small carbon enrichment is more likely to be attributed to the weak 
TDU in the primary star at such low metallicity (see Suda & Fujimoto, 2008). Thus we interpret 
that this star has evolved from a binary system with a primary of mass M >~ 3M Q , where the 
upper mass limit is set by the condition that dredged-up carbon can survive the hot bottom 
burning that are expected to work for massive AGB stars. It is to be noted that the stars of 
mass M > 2.5 M , even if born of primordial gas, start the TP- AGB evolution in the helium 
shell polluted by metals accreted on the surface since the surface pollutants have spread all over 
the hydrogen-rich envelope at the second dredge-up when the surface convection penetrating 
into the helium core (Suda et al., 2004). 

In conclusion, the peculiar abundance patterns of light elements Li through Al, observed 
from the three HMP/UMP stars can be explained in terms of the nucleosynthesis in AGB stars 
of population III. The three stars represent all three possible ways for the metal-free and EMP 
stars of low- and intermediate-mass to acquire the carbon-enrichment through helium shell- 
flashes, formulated by Fujimoto et al. (2000); Case II with He-FDDM only, Case II' with He- 
FDDM plus TDU and Case IV of TDU only. This gives a basis of the first-star interpretation 
of these stars, in which the enhancement of carbon and light elements is transferred from 
the erstwhile AGB companions and the other metals including the iron group elements are 
accreted from the interstellar gas, polluted with the supernova ejecta after their birth. It is 
true, however, that the pristine metals of such small level with the scaled-solar abundances 
relative to the observed iron abundances have little effects on the evolution and nucleosynthesis 
of primary stars and these stars can develop the observed enhancement of light elements through 
the binary mass transfer. Accordingly, we cannot exclude the possibility that these stars are 
the second-generation of stars, born out of gas already polluted by ejecta from the supernova 
explosion of the first generation of stars. One of the ways to discriminate these two possibilities 
of iron sources is the s-process nucleosynthesis as discussed in Suda et al. (2004). As seen 
from Figure 8, the neutron exposure well exceeds r cxp ~ 3 — 4 mb _1 (A( 13 C/ 12 C)/0.01) 1//2 , and 
hence, the heavy s-process elements can be synthesized even when initially the pristine metals 
as seeds are initially absent for sufficiently large amount of mixing. The pristine iron, if exists, 
also undergoes the s-process nucleosynthesis and is mostly converted into lead with such large 
exposure. In actuality, the enhancements of Sr and the upper limit of Ba are observed for 
HE1327-2326. We will discuss this problem in a subsequent paper (Nishimura et al. 2009 in 
preparation). 
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5. Conclusions and Discussion 

A distinctive feature of metal-free (Pop. Ill) and extremely metal-poor (EMP) stars of 
[Fe/H] < —2.5 is the occurrence of hydrogen mixing into the helium convective zone during 
helium shell flashes in an early phase of TP-AGB phase, which promotes the neutron-capture 
reactions with 13 C formed via the capture of mixed proton by 12 C as neutron source. We have 
studied the nucleosynthesis of light elements up to sulfur ( 35 S), attendant with this hydrogen 
mixing and revealed the characteristics of AGB nucleosynthesis for Pop. Ill and EMP stars. 
We demonstrate that they are the effective sites for the synthesis not only of carbon but also 
of oxygen and other light elements including the neutron-rich isotopes and odd elements. The 
resultant yields are compared with the abundance patterns observed for the EMP stars with 
carbon enhancement in the Galactic halo. The present results serve to identify the trace of AGB 
nucleosynthesis in Pop. Ill and EMP stars. It is also useful to single out the nucleosynthetic 
signatures of other mechanisms, in particular, of the first supernovae. 

Our main conclusions are summarized as follows; 
(1) Under the neutron-rich environment, the neutron and a-captures cooperate effectively to 
synthesize a variety of elements, starting from 12 C. The combination of (n,7) and (a,n) 
reactions produce 16 0, 20 Ne, and 24 Mg, in much more abundances than the (0,7) reactions. 
Neutrons, released via 13 C(«,n) 16 0, are spent mostly for 22 Ne production, and some for 26 Mg 
in the case of the strong flashes of large peak temperatures. The neutron capture reactions 
proceed with 22 Ne and 26 Mg as seeds to form heavier elements including odd nuclei of sodium, 
aluminum and phosphorus. The amounts of yields are determined essentially by the amount of 
mixed 13 C, with stronger dependence for heavier elements. The abundance of 16 increases in 
proportion to a square root of the amount of mixed 13 C until it reaches a significant fraction of 
12 C, while the abundance of 22 Ne is in proportion to the amount of mixed 13 C and attains at a 
half of mixed 13 C in number. There is a general tendency of yields decreasing more rapidly for 
heavier elements while local peaks form in the abundance distribution at the nuclei with the 
neutron capture cross section much smaller than those for the ambient nuclei as 22 Ne, 26 Mg 
and 31 P. 

(2) For the three most iron-poor stars, their peculiar abundance patterns from carbon 
through aluminum can be reproduced in terms of three possible modes of AGB nucleosynthesis 
of Population III stars born of primordial gas for the primary stars of different mass range: in a 
star of mass range M < 1.5 M & with helium shell flashes with hydrogen mixing and dredge-up 
by helium-flash driven deep mixing (He-FDDM) for HMP star of HE1327-2326 ([Fe/H] = -5.4): 
in a star of mass range ~ 1.5 < M\ < 3M with the additional nuclear products of helium shell 
flashes without hydrogen mixing and third dredge-up (TDU) following He-FDDM for the other 
HMP star of HE0107-5240 ([Fe/H] = -5.3): and in a star of mass in the range M > 3M with 
TDU alone for UMP star of HE 0557-4840 ([Fe/H] = -4.75). 
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(3) The oxygen enhancements, observed for carbon enhanced EMP (CEMP) stars, are 
well explained in terms of the nucleosynthesis in the helium shell flashes with hydrogen mixing 
and subsequent helium flash-driven deep mixing (He-FDDM) or even without hydrogen mixing 
and subsequent third dredge-up (TDU). This is the case in particular for all the CEMP-s stars 
with s-process enhancement [Ba(+Pb)/Fe] > 1 (CEMP-s stars), while some CEMP-nos stars 
without s-process enhancement [Ba/Fe] < 1 (CEMP-nos) exhibit the oxygen enhancement ex- 
ceeding those predicted from the AGB nucleosynthesis. As for the sodium, magnesium and 
aluminum, the observed abundances of CEMP-s stars are accommodated within the neutron- 
capture nucleosynthesis during helium shell flashes with hydrogen-mixing along with the con- 
tribution of pristine magnesium, while CEMP-nos lack the symptoms of the profound neutron 
captures. 

The agreement of the abundance characteristics for three HMP/UMP stars give the 
basis of our Pop. Ill binary interpretation of their origins, i.e., these stars emerge from the 
low-mass members of binary systems of the first generation stars with the different masses, 
suffering the surface pollution by accreting the wind matter from the primary AGB stars and 
also by accreting interstellar matter enriched with iron and other metals ejected by supernovae 
after their birth. The three stars represent the entire possible paths to carbon-enrichment in 
the binary scenario, as formulated in Fujimoto et al. (2000, see Suda et al. 2004, and Komiya 
et al. 2007); HE1327-2326 for case II, HE0107-5240 for case II', and HE0557-4840 for case IV. 
Our interpretation is consistent with the depletion of lithium in these stars. It is true that the 
pristine CNO elements, if they exist as much as the scaled-solar abundances with the observed 
iron abundances of [Fe/H] = —5.4 — 4.75, have nothing to do with the AGB nucleosynthesis 
as long as the light element synthesis is concerned. And yet it is worth pointing out that the 
degree of the surface iron pollution is in accordance with our assignment of the masses to their 
primary stars since the Bondi accretion rate increase in proportion to the second power of total 
mass. 

One of the merits of our scenario is that it is worked out within the current standard 
framework of stellar evolution in spherically symmetry; the stars of assigned masses must have 
brought about the characteristic abundances observed from these stars as long as they were 
formed. The prevalence of binarity is known for younger populations and we may take it for 
granted among the stars of smaller metallicity (Komiya et al., 2007). This forms a contrast with 
the supernova scenario, which assumes very artificial mixing-and-fallback model for supernova 
explosion, proposed by Umeda & Nomoto (2003). Recently, Tominaga et al. (2007); Tominaga 
(2009) attempt to interpret the mixing-and-fallback model as alleged jet-induced explosion, 
although it largely reduces the amount of carbon ejected and available for the formation of the 
next generation stars. For these models to work, there should be a very efficient mechanism to 
confine most of the carbon ejected by a supernova and form a star in a cloud of rather small 
masses (see, e.g., Suda et al., 2004). We also point out that the evolution of rotating stars 
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lacks the sound foundations since a reliable theory is yet to be properly established for the 
hydrodynamical instability and elemental mixing in rotating stars. 

In the present study, we treat the mixing process parameterized with the amount of 
mixed 13 C and the duration of mixing epoch. As seen in eq. (3), the mixing of 13 C, and hence, 
of hydrogen, is a small amount and at low rate not to cause the splitting of convective zone and 
He-FDDM. Such weak hydrogen mixing is likely to be missed in actual numerical computations 
since the entropy in the hydrogen burning shell starts as the helium flash convection keeps 
increasing after the peak of shell flashes (Fujimoto, 1977). In addition, the detailed modeling 
of hydrogen mixing will demand the proper treatment of the structure of convective motions 
and the entrainment of mixed hydrogen by the convection, which may be beyond the current 
framework of modeling of stellar structure and evolution. Rather, we may gain an insight into 
the proper treatment from the comparison with the observations of EMP stars. The same 
prescription is also applicable to the hot bottom burning. It is argued that the hot bottom 
burning grows efficient for lower metallicity, and some evolutionary calculations predict that 
carbon and oxygen are deeply processed by CNO cycles (e.g., Campbell & Lattanzio, 2008). 
On the other hand, it is true that there are only a few stars observed with the large N/C ratios 
(Johnson et al., 2007). Since the temperature at the bottom of surface convection depends 
sensitively on the treatment of convection, and hence, the constraints may be imposed from 
the comparison with the observations of EMP stars themselves. The effects of such internal 
mixing, whether they can be treated within the current standard framework of stellar evolution 
and beyond it, wait for future works, and have to be explored by recourse to comparisons with 
detailed abundances observed for individual stars. The present results may contribute to the 
theoretical effort in this line providing a basis for the understandings the characteristics of EMP 
stars. 
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Appendix 1. Neutron-Recycling Reactions and Oxygen Production in Metal-Poor 
AGB Stars 



Y 16 — ^12,0" 



exp{-a„ 7 ( ib 0)r} + F 16 ,oexp{-a n7 ( ib O)r}. (A2) 



If we focus on the allocation of neutrons produced by the (a,n) reactions, the progress 
of neutron-recycling reactions may be described in the form 

^ = Y 12 a n ,( 12 C) - F 16 a n7 ( 16 0) (Al) 

GST 

as a function of the neutron exposure, r. Here we may leave out 13 C since it takes a minor 
part as the destination of neutrons although the timescale is determined by the 13 C(a,n) 16 
reaction. Since 12 C is much more abundant than other nuclides in our case, we may well neglect 
the variation in Y12, and eq. (Al) yields an approximate solution analytically; 

^n 7 ( 12 C) r 

W 18 o) [ 

For r<l/ cr„ 7 ( 16 0), then, the production of 16 increases linearly with the neutron exposure as 
5Yiq = Yi 6 — Y"i 6i0 exp[— <T n7 ( 16 0)r] ~ Y"i2,o0Vy( 12 C)T, which is true mostly in our computations. 

The flow of nucleosynthesis bifurcates on the passage through 17 owing to the compe- 
tition between the neutron and a-captures, which may vary with the rate of hydrogen mixing 
and with the strength of shell flash. We denote the ratio between the neutron and a-capture 
rates of 17 by £ and the branching ratio of (0,7) and (a,n) reactions of 17 by /, i.e., 

C = Y 4 {(a an ( 17 0)v) + (a a ^ 17 0)v)}/Y n (a na C 7 0)v), (A3) 
/ = (cr a ,( 17 0)v}/ ((a an ( 17 0)v) + (a a ,( 17 0)v)) . (A4) 
Then, we may approximate the further progress of nucleosynthesis up to 22 Ne as 



Y w a ni ( lb O) - Y 17 a na ( u O)(l + 0, (A5) 



dYn _ 
dr 

^ = Fi7^( 17 0), (A6) 
dr 

Y l7 a na ( 17 0)^ ~ f) ~ >W( 20 Ne), (A7) 



dr 
dY 

~ Y l7 a na { l7 0)if + F 20 a n7 ( 20 Ne) - F 21 a n7 ( 21 Ne), (A8) 



dr 

>2ia n7 ( 21 Ne) (A9) 



dY 22 , 2 i 



dr 

where we neglect the absorption of neutron by 14 C because of smaller cross section. For the 
weak shell flashes of longer timescales, 14 C may decay into 14 N that has much larger neu- 
tron capture cross section. This makes little difference, however, since 14 N reacts mainly 
through (n,p) channel and the released protons are converted again into neutrons through 
the 12 C(p,7) 13 N(e+z/) 13 C(a,n) 16 reactions. 

If we regard £ and / as constant, neglecting the variations of the reaction rates with 
the temperature, the above system of equations yields an analytical solution. Figure 12 shows 
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the evolution of number abundances of ' O, C and 20 ' 21 ' 22 Ne for ( = 1 as a function of 
neutron exposure; this reproduces the characteristics observed from the numerical computations 
in Figure 5 very well. The abundance of 16 increases linearly with the neutron exposure, as 
stated above, and approach an upper bound of Y\^jY\ 2 = a n7 ( 12 C)/cr n7 ( 16 0) for r > l/o" n7 ( 16 0). 
In comparison, Yn and Y\± remain in much smaller abundances. They increase with r in 
proportion to higher powers as 

Y 17 ~ <r n7 ( 16 O)[Yi 6 , T + (l/2)F 12i0 a n7 ( 12 C)r 2 ] (A10) 
Y u ~ a na ( 17 O)[(l/2)F 16 , a n7 ( 16 O)r 2 + (l/6)F 12 , a n7 ( 12 C)a n7 ( 16 O)r 3 ] (All) 

for small r. For r > l/cr na ( 17 0)(l + £), the former tends to reach a steady-state value, Yyj ~ 
^i2C r n7( 12 C)/cr riQ ,( 17 0)(l + £). On the other hand, F14 keeps increasing simply because of the 
neglect of a-capture of 14 C. In actuality, 14 C as well as the daughter nuclei 14 N eventually burn 
into 18 and into 22 Ne since the a-captures of these nuclei and of 18 are much faster than the 
3a reactions (see Figure 4). 

The production of 20 Ne and 21 Ne is in proportion to £. Between them, 20 Ne is much 
more abundant than 21 Ne because of larger branching ratio (1 — /)// and smaller neutron 
capture cross section of 20 Ne. For large neutron exposure, their abundance ratio approaches a 
steady-state value of Y" 20 (l — /O^Vy^Ne) ~ Y^ov^^Ne), and eventually, 22 Ne overweighs them. 
Note that 22 Ne plotted in this figure is only those formed solely through neutron captures while 
much more is formed via the a-captures of 14 C. 

In the recycling reactions, neutrons are spent mostly to form neutron-rich isotopes of 
oxygen and neon so that the progress of nucleosynthesis is limited by the number of available 
neutrons, and hence, by the amount of mixed 13 C. Accordingly, the termination is given by the 
relation 

AF 13 , m ix = Y 17 + Y 21 + 2 x (Y 14 + Y 22 ). (A12) 

More neutron may be released from 17 and neon isotopes. The former is effective when the 
a-capture overweighs the neutron capture and the latter are limited to the shell flashes of high 
temperatures (T > 3 x 10 8 K). 

Figure 13 shows the final abundances of 16 and neon isotopes, 20 Ne and 22 Ne, where 
all 14 C are assumed to be converted into 22 Ne by the a-captures. The oxygen production 
increases with the mixed amount of 13 C and reaches as large as a half of the upper bound, 
Y 16 /Y 12 = CT n , 7 ( 12 C)/(x n7 ( 16 0) ~ 0.405 for A( 13 C/ 12 C) > 0.01. The amount of oxygen produced 
is limited by the leak of neutrons from the recycling reactions as they are captured by 16 0. 
It is followed by the capture of another neutron, while a fraction, £(1 — /)/(l + £) of neutrons 
are returned by the 17 O(a,n) 20 Ne reaction. Thus, the condition that all the neutrons leave the 
recycling reactions may be written in the form; 

AF 13 , mix = J a„ 7 ( 16 0)F 16 |2--^y(l-/)|dr. (A13) 
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This integral can be carried out with use of eq. (A2), and in particular, for r <C l/er n7 ( 16 0), 
the condition reduces to; 



As the amount of mixed 13 C increases, then, the second term in the right-hand side member 
overweighs the first term, and the abundance ratio of O to C increases in proportion to a 
square root of the ratio, A( 13 C/ 12 C), between the mixed amount, AYi3 im ; x , of 13 C to the 
abundance Y\2 of 12 C. In this figure, we also plot the average number, N cyc , of recycling reactions 
12 C(n, 7) 13 C(a, n) 16 that each neutron experiences, estimated by = 5Yi e / AY 13 mix . It 
increases with A Yj 3 ) TT1 j x from N cyc ~ Yi 20 a ni ( 12 C) /Y 16t0 a ni ( 16 O) to reach a maximum, which is 
larger by a factor, (1 + £)/(l + £/), for larger £ owing to the contribution of (a,n) reactions, 
and then, turns to decrease inversely proportional to the square root of A ( 13 C/ 12 C) because of 
the increase of 16 O. 

The abundance of 20 Ne depends on the ratio, £, of neutron to a-capture rates of 17 
while 22 Ne is in proportion to mixed amount of 13 C and Y22 — 0.5 Alia mix- As a consequence, 
20 Ne can overweigh 22 Ne for large £, i.e., for a strong flash of higher temperature, while the 
neutron capture grows important to convert 20 Ne into 22 Ne for a larger exposure of r > 8.4, 
and hence, for large mixing amount of A( 13 C/ 12 C) > 0.04. 

Appendix 2. Surface Sodium Enhancement during Third Dredge-up 

When the third dredge-up occurs, the carbon abundance in the envelope may greatly 
increase. A part of the dredged-up carbon is burnt in the hydrogen burning shell to be 14 N, 
and then, capture a particles in the helium flash convective zone to be 22 Ne, which is brought 
back to the envelope by the succeeding third dredge-up. Furthermore, 22 Ne, thus enriched, is 
burnt in the hydrogen burning shell and some is converted to 23 Na, and again, brought back to 
the envelope by the third dredge-up. 

If we denote the carbon abundance in the envelope after the i-th dredge-up by 
Y12, the surface abundance of 22 Ne is increased by the z-th dredge-up as much as 5Y 22 = 
(X12 /^.He.conv)^^ while the surface carbon abundances is increased by SY^ 2 (where Yi 2) H e ,conv 
is the carbon abundance in the helium flash convection). The increment of sodium abundance in 
the envelope is given by SY 23 = giY^ 1 /Y 12 ^ conv )5Y^ 2 , where g is the fraction of 22 Ne, converted 
to 23 Na by the NeNa chain reactions. In the limit of small variations, we may approximate the 
increases in the abundances by the differential equations, which yields the analytic solutions; 




(A14) 



Y 22 — Y 22 fl + (Y l2 — Y 120 )/2Y 12 



(A15) 



_ \ 1 12,Hc,conv 12,He,conv/ u 1 12,He,conv 

Here the subscript denotes the abundance at the start of third dredge-up. 





(A16) 
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Since the contribution of pristine metals is negligible in the EMP stars, the sodium 
enrichment of AGB stars may be related to the carbon abundance as; 

[Na/C] = -2.2 + log^ + 2[C/H] - 21og(F 12 , H c,co„v/0.017). (A17) 

The value of g decreases rapidly with the temperatures in the range relevant to the hydrogen 
shell-burning; it varies from g ~ 1 for T < 5 x 10 7 K, to g~ 0.1 at T ~ 6 x 10 7 K, and to g ~ 0.01 
at T > 7 x 10 7 K (Angulo et al., 1999), although the reaction rates of NeNa cycles are subject 
to fairly large uncertainties (some by more than a order of magnitude El Eid et al., 2004). 

In the massive AGB stars, carbon is burned in the envelope convection (HBB) to be 
converted into nitrogen, which enlarges the [Na/C] ratio; this is argued to be the case for stars 
with EMP stars without the carbon enhancement such as "mixed stars" (Spite et al., 2005), 
where [Na/C] varies from —0.5 to 1.0. Furthermore, for still more massive AGB stars, it is 
possible to burn neon by HBB, which may raise the sodium enrichment up to; 

[Na/C(+N)] ~ -0.07 + log<? + [C(+N)/H] - log(F 12 , Hc ,conv/0.017). ( A18) 

If the temperature during the shell flash reaches above ~3x 10 8 K, on the other hand, 
most of 22 Ne is burnt into magnesium in the helium flash convective zone. For such massive 
stars as enter into TP-AGB phase with the sufficiently large core, we expect the magnesium 
enrichment as large as 

[Mg/C] = -1.4 + [C/H] -log(y 12) He,conv/0.017). (A19) 

In this case, massive isotopes of 26 Mg and 25 Mg are enhanced. 

The sodium enrichment is observed at [Na/C] = —1.7 for HE 1327-2326 and at much 
smaller [Na/C] = -3.2 for HE 0107-5240. In the above equation of eq. (A17), [C/H] is the 
carbon enrichment in the primary stars in our binary scenario and can be larger than currently 
observed for the low-mass secondary member, polluted via wind accretion. If [C/H] ~ — 0.5, 
the sodium enrichment of HE0107-5240 can be consistent with that by the third dredge-up 
alone for g = 0.1 — 0.01. It is difficult, however, to reproduce the observed sodium enrichment 
for HE 1327-2326 without hydrogen mixing and He-FDDM. However, the sodium enrichment 
of HE 0107-5240 can be consistent with the enrichment in eq. (A17) with [C/H] ~ — 0.5, 
although in this case seek we may not expect the enrichment of aluminum. The rather low 
enrichments of magnesium observed for HE0107-5240 and for HE 1327-2326 ([Mg/C] = —3.85 
and —2.5, respectively) are indicative of low-mass nature of their erstwhile primary stars. For 
HE0557-4840, [Na/C] = —1.81 and [Mg/C] = —1.4, and yet, the sodium and magnesium are not 
enriched in comparison with iron ([Na/Fe] = —0.16 and [Mg/Fe] = 0.25) as well as aluminum 
at [Al/Fe] = -0.61 (Norris et al., 2007). It is also noted that HE0557-4840 is distinguishable 
from HE0107-5240 and HE1327-2326 in the degree of carbon enrichment, which is indicative of 
the difference in their erstwhile primary AGB stars and in the efficiency of third dredge-up. 
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Table 1. Model Parameters 
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Y 

16, p 


1 
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2.046e-2 


20.001 
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4.33 


8.758e-3 


1.141e-5 


1.205e-2 


3.596e-5 


2 
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1.631e-2 


19.821 


8.348 
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1.761e-2 


4.449e-4 


2.137e-2 


4.973e-4 
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5 
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1.307e-2 


20.881 


8.521 


9.13 


1.716e-2 


2.506e-4 


2.193e-2 


2.575e-4 



a : mass and radius at the bottom of helium burning shell 
b : proper pressure of the helium shell-flash 
c : maximum temperature reached during helium shell flashes 
d : maximum helium burning rate at the peak of shell flashes 
e : mole abundances of 12 C and 16 O at the onset of shell flashes 
f : mole abundances of 12 C and 16 O at the peak of shell flashes 
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Fig. 1. A schematic drawing of the progress of mixing during helium shell flashes in extremely metal-poor 
stars. The convection driven by a helium flash extends outwards beyond the bottom of hydrogen-containing 
layer; mixed hydrogen is carried down into the convective zone to be captured by 12 C and to form C 
after /3-decay. 13 C is further mixed down to the bottom of helium convective zone and release neutrons to 
initiate neutron-capture reactions and synthesize s-process nucleosynthesis (left panel). When the helium 
shell flash grows strong and mixes hydrogen at a sufficiently large rate, eventually, the burning of mixed 
hydrogen splits the flash convective zone into two, i.e., an upper one driven by the hydrogen shell flash and a 
lower one driven by the helium shell flash. During the decay phase of the shell flash, the surface convection 
penetrates into the shells involved in the upper convective zone to cause helium-flash driven deep mixing 
(He-FDDM) and carry out nuclear products to the surface (middle panel). For stars of M > 1.5 Mq, third 
dredge-up events (TDU) occur during the following shell flashes. Carbon and oxygen, newly synthesized 
by 3a reaction and a-capture of 12 C during the helium shell flashes without the hydrogen mixing, are 
dredged up to the surface along with the remainder of nuclear products, formed by the former neutron 
capture burning and stored in the helium zone (right panel). For the detailed description about the model, 
see the following references (Hollowell et al., 1990; Fujimoto et al., 2000; Iwamoto et al., 2004; Suda et al., 
2004; Komiya et al., 2007). 
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Fig. 2. The structural changes during the helium shell flash, calculated with use of the analytical solution 
of helium flash in finite amplitude (see e.g., Fujimoto et al., 1999); top and bottom panels show the time 
variations of temperature and density at the bottom of flash convcctivc zone. 
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Fig. 3. Flowchart of nucleosynthesis during the convective J C burning. Right arrows denote the neutron 
captures to which the neutron capture cross sections are attached in unit of mb at thermal energy kT = 30 
kcV, taken from Bao et al. (2000). Upward arrows denote the (a, 7) reactions and upward and downward 
diagonal arrows denote the (a,n) and (n,a) reactions, respectively. Thick arrows represent the main paths 
of nucleosynthesis. 
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log T (k) 

Fig. 4. The rates of a-capture reactions d\ogYi/dt = Y^pN a{&v) a -y and Y^pNAivvjan, relative to the 
3a reaction rate (dlogl^/ 'dt = 3Y^p 2 N^(av) aaj /6), as a function of temperature. Rates are taken from 
Angulo et al. (1999) except for the rates of 14 C(a,7) 18 0, 16 O(a,7) 20 Nc, 17 0(a, 7 ) 21 Nc, 21 Nc(a,7) 26 Mg, 
24 Mg(a,7) 28 Si and 25 Mg(a,n) 28 Si taken from Caughlan & Fowler (1988). We adopt 3a reaction rate from 
Caughlan & Fowler (1988) for the a-capture reaction rates, enumerated above. The density and helium 
abundance are set at p = 10 3 g cm -3 and Y4 = 0.2, respectively. 
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Fig. 5. The progress of nucleosynthesis, induced by hydrogen mixing, in the helium-flash convection. 
The elemental abundances are plotted against time (left panel) and against the neutron exposure, r (right 
panel), respectively; the neutron abundance, X n , is multiplied by 10 13 . The results are taken from Model 
4 with the amount of mixed 13 C, AXi3. m ; x = 0.01Xi2(t pca k) in total, (where Xi2(i pca k) is the abundance of 
12 C in the helium convective zone at the peak of flash, t = i p0 ak), or A( 13 C/ 12 C) = 0.01, which is supplied 
at a constant rate for an interval of At m ; x = 10 9 s starting from the peak of shell flash. 
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Fig. 6. The number abundances of nuclear products, normalized by carbon abundance, are plotted as a 
function of mass number, A, resulting from the nucleosynthesis in the helium shell flash convective zone. 
The results are shown for Model 1 with Ai mix = 10 12 s (top panel) and Model 5 with A< mix = 10 9 s 
(bottom), respectively. Each line denotes the yield for the different amount of mixed 13 C, A( 13 C/ 12 C) as 
indicated in top-right corner of each panel. 
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Fig. 7. The same as in Figure 6 but for the different durations, At m ; x , of mixing epoch with the same 
amount of mixed 13 C, A( 13 C/ 12 C) = 0.01 for Models 1 (top panel) and 5 (bottom panel): numerals in the 
top-right corner of each panel denote the durations of mixing epochs in units of second. 
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Fig. 8. The neutron exposure as a function of time for Model 1 (top panel) and Model 5 (bottom panel) 
with the different amounts and duration of 13 C mixing. The amounts of mixed 13 C, A( 13 C/ 12 C), are 
denoted by different lines as indicated in the top-left corner. The durations of mixing Ai m j X are read from 
a break in each curve. 
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Fig. 9. Comparison of the elemental abundances of yields of AGB nucleosynthesis resultant from hydro- 
gen mixing into the helium flash convective zone with the observed abundance patterns of light elements for 
three most iron-deficient stars; the computational results are plotted, normalized with 12 C abundance of 
HE1327-2326 for Model 1 with the different amounts of mixed 13 C, A( 13 C/ 12 C) = 0.0001, 0.001, 0.01, and 
0.02, 0.03, and 0.1 and all with At m i X = 10 12 s, by solid lines from bottom to top, respectively. The upper 
dotted line is the result for Model 3 of higher peak temperature with A( 13 C/ 12 C) = 0.02 and At m ; x = 10 10 
s. The lower dotted line is obtained by multiplying the above result by 1/20 under the assumption that 
the elements are diluted by a factor of 20 relative to the carbon abundance which increases due to the 
third dredge- up subsequent upon He-FDDM(see Suda et al., 2004). The observed abundances and upper 
limits are taken from Christlieb et al. (2004) for HE 0107-5240, from Aoki et al. (2006) for HE1327-2326, 
and from Norris et al. (2007) for HE 0557-4840. The oxygen abundance is taken from Bessell et al. 
(2004) for HE0107-5240, and the middle of three values (see Figure 10) taken from Frebel et al. (2006) for 
HE1327-2326, while only upper limit from Norris ct al. (2007) for HE 0557-4840. 
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Fig. 10. Comparison of the oxygen abundances, theoretically predicted from the nucleosynthesis in AGB 
stars, with the observations for HMP/UMP and EMP stars, where the ratios of the oxygen to carbon 
abundances are shown against the carbon enhancements; thick solid line denotes the upper limit of the 
[O/C] values, reached by hydrogen mixing and He-FDDM. Two broken lines bound the range that can 
be reproduced by TDU without hydrogen mixing. Symbols denote the observed oxygen abundances for 
stars of [Fc/H] < —2.3, taken from Stellar Abundance for Galactic Archeology Database (SAGA Database; 
Suda et al., 2008); open circles and square denote HE1327-2326 and HE0107-5240, respectively and di- 
amonds show the upper limits of oxygen abundances derived under two different effective temperatures 
for HE0557-4840 (Norris et al., 2007): filled triangles and open triangles denote CEMP-s and CEMP-nos, 
respectively, while filled and open pentagons show those based on the medium-resolution spectroscopy 
by Beers et al. (2007): and crosses denote EMP stars without carbon enhancement ([C/Fe] < 0.5). For 
HE1327-2326, plotted are three oxygen abundances, an upper bound derived for OI triplet line A777 nm 
(middle), an abundance by ID LTE analysis of UV OH features (the highest one), and one with 3D correc- 
tion (the lowest one) (Frebel et al., 2006). We show two data by ID LTE analysis of UV spectra and with 
3D corrections included for HE1300-0157 (right and left ones, respectively Frebel et al., 2007), the upper 
limit derived from [OI] A6300 as well as the abundance from OI triplet for LP 625-44, (lower upper ones, 
respectively Aoki et al., 2002), and LTE oxygen abundances from OI triplet line A777 nm and that with 
NLTE correction for CS31080-095 (Sivarani et al., 2006). For CEMP stars with the nitrogen enrichment 
larger than the carbon enrichment, the abundances where nitrogen is added to carbon are indicated by 
open triangles linked by dotted line. Two oblique thin solid lines denote the loci of [O/Fc = 0.5] and 1.0, 
and thin broken lines delineate the ranges for a mixture of the yields of AGB nucleosynthesis and the 
envelope matter of the pristine oxygen abundance [O/C] = 0.5 and 1.0. 
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Fig. 11. Abundance ratios among sodium, magnesium and aluminum, predicted from the nucleosynthesis 
of neutron capture reactions in the helium flash convection attendant upon hydrogen mixing and the 
comparisons with the observations for HMP/UMP and EMP stars on the [Mg/Na] and[Al/Mg] diagram. 
Thick solid lines with arrow depict the evolutionary paths for some cases and small circles denote the 
end products of numerical computations. Labels (Ml , M4, and M5) attached to solid lines indicate the 
results for Model 1 with A( 13 C/ 12 C) = 0.1 and A< raix = 10 12 s, for Model 4 with A( 13 C/ 12 C) = 0.1 and 
At m ix = 10 10 s and for two Models 5 with Ai m j X = 10 6 s (left) and Ai m ; x = 10 9 s (middle with upward arrow) 
and both with A( 13 C/ 12 C) = 0.1, respectively, and broken and dash-dotted lines connect the sequences 
of different At m ; x with fixed A( 13 C/ 12 C) = 0.01 for Models 3 (lower) and 5 (upper), respectively. The 
observed abundances are collected with use of SAGA Database and denoted by the same symbols as in 
Figure 10. Two horizontal and vertical dotted lines indicate the upper limits to the abundance ratios, 
[Mg/Na] and [Al/Mg], that can be reached solely by the neutron capture reactions. 
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Fig. 12. Progress of nucleosynthesis in the helium flash convection with the neutron capture reactions, 
plotted against the neutron exposure r. Variations of abundances are computed with use of the system of 
approximate equations given in Appendix 1 by assuming that the ratio between the neutron and a-capturc 
rates of 17 equals to unity, i.e., £ — 1 [sec eq. A3]; thin solid labeled 5 16 denotes the production of 16 0, 
SYia = Yig — Fi6,o cx p[— ct„ 7 ( 16 0)t]. Note that the a-captures of 14 C is neglected, and if they are included, 
14 C will be converted into 22 Nc. 
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Fig. 13. Production of oxygen and neon by the doubly neutron-recycling reactions as function of the 
mixed amount, A( 13 C/ 12 C), of mixed 13 C relative to 12 C [= AXi3 im i x /Xi2(t pe ak)] for the different ratio, 
£, between the a and neutron capture reactions of 17 0. Broken, solid, and dash-dotted lines denote 
the results for £ = 0.1, 1, and 10, respectively. Dotted lines on the top denote the number, A^yc, of 
neutron recycling reactions 12 C(n,7) 13 C(a,n) 16 that each neutron experiences, and hence, the ratio of 
the produced oxygen to the mixed 13 C (<5 16 0/AYi3 jm j x ) for £ = 0.1, 1, 10, from bottom to top. 
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